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(54) Precede de localisation spatio-temporelle, de radio-mobiles en milieu urbain 

(57) Le procede consiste : 



a calculer (1) la reponse impulsionnelle du canal de 
propagation entre le radio-mobile et la station emet- 
trice/receptrice, 

a estimer (2) a partir de la reponse impulsionnelle 
les couples (angles, retard) du trajet de propagation 
direct et des multitrajets resultats de reflexions sur 
des obstacles situ6s dans le champ de propagation 
des ondes entre la station de base et le radio- 
mobile. 

a calculer (2) pour chaque trajet les attenuations otj, 
et a calculer (3) la distance d s6parant la station 
emettrice/receptrice du radio-mobile en conside- 
rant seulement les retards et les attenuations du 
trajet direct et d'un multitrajet quelconque. 

Applications : Localisation et suivi de radio-mobiles. 
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Description 

[0001] La presente invention concerne un procede de localisation de radio-mobiles en milieu urbain. 
[0002] Elle s'applique notamment a la localisation et au suivi de radio-mobiles dans un reseau GSM pour anticiper 
leur passage d'une cellule a la suivante et ameiiorer la reception des stations de base par la connaissance de la posi- 
tion des radio-mobiles emetteurs. 

[0003] Dans des configurations simples de propagation, la localisation de radio-mobiles est traditionnellement obte- 
nue a I'aide de radio-goniometres disposes en des lieux determines d'un territoire. La position d'un radio-mobile peut 
alors etre determine par triangulation en reportant sur une carte geographtque les lieux ou sont situes les radiogonio- 
metres ainsi que les directions d'arriv6e sur ceux ci des ondes radioeiectriques emises par le radio-mobile. 
[0004] Cependant, cette methode qui repose sur I'exploitation des seuls chemins de propagation directs apparaTt dif- 
ticilement applicable dans le cas ou les ondes empruntent plusieurs trajets et est egalement lourde de mise en oeuvre 
car elle necessite plusieurs stations de reception. Ce type de propagation est notamment celui qui est rencontre en 
milieu urbain dans les r6seaux de communication radio-mobiles cellulaires, ou le spectre des communications est situe 
dans une gamme d'ultra hautes frequences s'etendant typiquement de 30 MHz a 3 GHz pour.le reseau GSM. Pour 
cette gamme de frequences la propagation des ondes radios ectriques entre un emetteur et un recepteur peut 
empreinter des trajets multiples dus aux reflexions et aux diffractions sur des obstacles interposes sur le chemin de pro- 
pagation direct entre remetteur et le recepteur. La gamme des frequences utilisees est aussi la cause de trajets diffus, 
aussi bien temporellement que spatialement, qui degradent les performances des liaisons. 
[0005] Le but de ('invention est de pallier les inconvenients precites. 

[0006] A cet effet. invention a pour objet un procede de localisation spatio-temporelle d'un radio-mobile relativement 
a une station de base emettrice/r6ceptrice caracterise en ce qu'il consiste : 

- a calculer la reponse impulsionnelle du canal de propagation entre le radio-mobile et la station emettrice/recep- 
trice, 

- a estimer a partir de la reponse impulsionnelle les couples (angles, retard) du trajet de propagation direct et des 
multitrajets resultats de reflexions sur des obstacles situes dans le champ de propagation des ondes entre la sta- 
tion de base et le radio-mobile. 

a calculer I'attenuation aj de I'onde propagee par chaque trajet, 

- et a calculer la distance d separant la station emettrice/receptrice du radio-mobile en considerant seulement les 
retards et les attenuations du trajet direct et d'un trajet non direct quelconque. 

[0007] D'autres caracteristiques et avantages de ('invention apparaitront a I'aide de la description qui suit faite en 
regard des dessins annexes qui represented : 

la figure 1 le procede selon I'invention sous la forme d'un organigramme, 
la figure 2, un modeie de reception monocapteur, 

- la figure a, une illustration d'une reception multicapteur par une antenne circulate, 

les figures 4 et 5, une illustration d'une transmission a deux bu trois trajets en milieu urbain. 

[0008] Le procede selon I'invention se deroule de la fagon representee a la figure 1 suivant trois etapes principales. 
Dans une premiere etape, representee en 1 a la figure 1 , il consiste a calculer la r6ponse impulsionnelle h(t) du canal 
de propagation pour estimer a retape 2 les angles 6 m et les retards x m des differents trajets ainsi que les attenuations 
otp des ondes propag6es sur chacun des trajets et enfin a calculer a I'etape 3, a partir des retards x m calcuies et de la 
reponse impulsionnelle h(t) du canal, la distance d separant un radio-mobile de la station demission/reception. 
[0009] Le calcul de la reponse impulsionnelle du canal qui est effectue a l'6tape 1 a lieu sur le signal obtenu a la sortie 
d'une antenne multi-capteurs. 

[001 0] En presence de plusieurs trajets M d'emission bande etroite de radio-mobiles, le signal obtenu a la sortie d'un 
capteur n de I'antenne multi-capteurs est 6gal a la somme des M signaux s m (t) 6mis dans les directions 6 m retard6s de 
x m et multiplies par la reponse aj(6 m ) du capteur i dans les directions e m a laquelle s'ajoute un bruit bj(t). Dans I'hypo- 
these d'une transmission num6rique a bande etroite de frequence, une representation de la sortie vectorielle de 
rantenne compos6e de N capteurs s'obtierrt en concatenant la sortie x t (\) de chaque capteur. chaque sortie x t {\) veri- 
f iant la relation : 

M 

x,(t)= ^s^-a^ej+b.O) (1) 

m=1 
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[0011] En introduisant un vecteur directionne! a(6 m ) tel que 



^ml^a^ej a.tej] 1 (2) 



et en concatenant la sortie Xj(t) de chaque capteur, on obtient un vecteur x(t) tel que : 

M 



x(t) = [x^t), Xj (t), x N (t)] 1 = £a(6 m )s m (t) + b(t) ( 3 ) 



[0012] En rassemblant les vecteurs directionnels a(e m ) dans une matrice 

A = [a(e 1 ) i a(e 2 ).-..a(e N )] 

la relation precedente s'ecrit encore : 

x(t)=As(t)+b(t) (4) 



avec s(t) = [ Sl (t).s 2 (t). • •s M (t)] , 1 
ou St est le signal du 1 er trajet, 
s 2 est le signal du 2 ieme trajet, 

s M - est le signal du M ieme trajet. 

[0013] Suivant le modele de propagation qui est represents schematiquement a la figure 2, le signal x(t) qui est 
obtenu a la sortie d'un capteur de reception est un signal de la forme : 

*W= £« m .e , * m .a(8 m ).s(t-T m ) + b(t) (5 ) 

m=1 

ou a m designe (attenuation du m i6me trajet 
i m est le retard de groupe du m teme trajet 
M est le n ombre de trajets 
et <t> m est le dSphasage du m i6r ™ trajet. 

soit encore sous la forme d'une convolution du signal 6mis s(t) avec la r6ponse impulsionnelle h(t) du canal : 

soft x n (t)«h n (t)* sffl + b n (t) (6) 

pu b n (t) represente Tenveloppe complexe du bruit additif. 
[0014] Dans le cas d'une reception sur I'antenne composee de multi-capteurs, comme cela est represents k la figure 
3, la reponse impulsionnelle du canal est modifiee par la reponse de I'antenne suivant les directions d'arrivSe. Cette 
reponse est de ia forme : 

M 

h n (t)= £a n (e nvYm )a m e i H(t-T m ) (7) 

m=1 
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en supposant le canal ideal, c'est-a-dire de bande inf inie. 
[0015] En mode 6chantillonne, la r6ponse du canal de bande B prend la forme : 



h n [k] = F e £ a n (6 J a m e sin c((k - F e x J B) 



(8) 



m=1 



avec k € Z, Tespace des nombres entiers. 



10 



et sinc((k- F e i m ) B) = 



sin(7i(k-Fe-c m )B) 
n(k-F e T m )B 



(9) 



15 F e est la fr6quence d'echantillonnage 

B est la bande en fr6quence r6duite du signal. 

[001 6] Cette reponse pert §tre estimSe par un f iltre causal, a reponse impulsionnelle f inie de longueur L et de coef- 
ficients tels que : 



20 



h n [k] = £ a n (6 m )ct m e '* m sin c((k - F e x m ) B) 



(10) 



25 



pour ke {0 ... L} 

et h n [k] =0pourk <0ouk )L 



30 



35 



Le vecteur h n correspondant s'6crit h n = I n n(°)' * * * ■ h n( L )] 1 



(11) 



[0017] Suivant un premier mode de realisation de ('invention, Testimation des angles e m et des retards x m de l'6tape 
2, est execute de facon disjointe, en utilisant la relation matricielle suK/ante : 

H = [h 1 ,h 2l ...,h N ] = M.aA t (12) 



avec : 



40 



45 



M = 



sincC-x^eB) sinc(-t 2 F e B) 
sinc((l - ?^Fe)V) sinc((l - x 2 F e)B) 



sinc(-T M F e B) 
sinc((l-i M F e )B) 



sinc((L - xf e ) B) sinc((L - x 2 F e ) B) sinc((L - x M F e ) B) 



(13) 



50 



55 



4 



EP0 926 510 A1 



£1 = 



<*1 
0 



0 

a 2 



: '■. '•. 0 
0 ••• 0 a N 

A = [a 1 (9 1 )...a N (9 N )] 



(14) 



(15) 



[001 8] (.'estimation des parametres temporels t m s'obtient a partir de la relation (1 2) en estimant un projecteur 

W 

sur lespace temporel HH* et en appliquant soit la methode connue sous labreviation anglo-saxonne MUSIC, 
de ( (Multi Signal Classification )>en cherchant les maxima de la quantite 

P t= _! (16) 



OU 



m(T) = [sinc^-xFeB) .... sinc((L - tF^B)] 1 



est une colonne de la matrice M(13), soit la methode deterministe connue sous le nom de < ( maximum de vraisem- 
blance )) en utilisant un critere de la forme 

J(t 1( t m )= t r ((l L+1 -M (M^-^tjHHt) (17) 



[001 9] Une description de la methode MUSIC peut §tre trouvee dans Particle de M. Ralph O. Schmitt ayant pour titre 
<< Multiple source DF signal processing : An experimental system >>publi6e dans IEEE transaction on antennas and 
propagation AP34 N° 3 March 1986. 

[0020] La determination des M couples (6 m , t m ) permet de localiser les radio-mobiles. Dans le cas ou, comme cela 
est repr§sente a la figure 4, seulement deux trajets, un trajet direct et un trajet comportant une reflexion sont seiection- 
nes, la determination des couples a lieu en considerant quatre parametres e 1( t 1( e 2 , t 2 qui se ramenent a la connais- 
sance de 6^ , e 2 et des retards , x 2 . En appelant d la distance du trajet direct qui sSpare remetteur radio-mobile 4 d'un 
rdcepteur 5 et d n , d 2 les distances du trajet indirect s6pararrt respectivement le recepteur 5 d'un obstacle 6 sur lequel 
I'onde indirecte est retiechie et remetteur 4 de I'obstacle 6, la distance d est fournie par la relation : 

d = i (18) 

ou a, et ct2 designent les attenuations du premier et du deuxieme trajet et \i la perte due a la reflexion du second 
trajet Les attenuations a, et sont obtenues a partir des couples (9 1f x, ; 6 2 , x 2 ) et de la relation matricielle 
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is 



25 



30 



CI = M H A 



a -i 0 
0 a 2 



(19) 



ou otj correspond a I'attenuation du i 6me trajet. 
[0021] Comme la perte de puissance dans une propagation libre se fait en 



_1_ 
r 2 



ou r est la distance parcourue, a-, et a 2 sont Ites par la relation : 



a, = (d^) 2 (20) 
a 2 d 2 

(T 2 -T 1 ) a .C 2 +2d(T 2 -C 1 )C 

avec d 1 = 2[(x 2 -x 1 )c^Kj-dcos(e 2 -e 1 )] 



et d 2 = ^d 2 +di -2dd,cos(e 2 - (22) 



[0022] S'il y a plus de deux trajets. la figure 5 montre que Information apportee par tout trajet supplemental se 
resume a deux parametres, et les nouvelles inconnues sont egalement au nombre de deux. En consequence, une 
detection de plus de deux trajets n'a pas d'inter&t dans le sens ou elle ne supprime pas les inconnues. II suffit dans ces 
conditions de limiter les traitements en considerant seulement le trajet direct et un trajet autre quelconque. Connaissant 
les angles d'arrivee des differents trajets, il est possible, grace & une carte topographique detailiee en trois dimensions 
de trouver les points de reflexion les plus probables et ainsi de remonter a la source demission. 
35 [0023] D'autres modes d'estimation des parametres sont egalement possibles. 

[0024] L'estimation disjointe des parametres d6crite pr^cedemment a I'avantage de reduire I'espace de recherche. 
Cependant. dans cette methode , les angles et les retards ne sont pas associ6s si bien qu'elie souffre d'un manque de 
resolution pour estimer les retards proches. 

[0025] D'autres modes de mise en oeuvre de rinvention permettent d'associer les parametres en estimant la matrice 
40 des attenuations il avec les difterentes combinaisons (e m , x m ) possibles en cherchant les couples rendant la matrice 
a la plus diagonale possible. 

[0026] Selon un deuxieme mode de mise en oeuvre de Hnvention, le precede op6re par estimations successives en 
estimant successivement les angles 0 et les retards x tout en les associant. Pour cela, deux m£thodes peuvent etre uti- 
lisees. La premiere consiste a estimer les parametres temporels x m connaissant les angles 6 m . 
45 [0027] En estimant des angles 6 m par la methode MUSIC ou la methode du maximum de vraisemblance l'estimation 
des retards x m est obtenue en estimant la matrice M par la relation matricielle 



M = HA * (A 1 A) 1 (23) 



50 



[0028] Chaque colonne de la matrice M forme un projecteur permettant d trouver le retard correspondant x m en 
appliquant lalgorithme MUSIC pour une seule source demission. 
55 [0029] Le critere utilise est de la forme : 



6 
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10 



15 



mix^Pi m(T) 

P(t) = ; ! (26) 

m r (x) rn(T) 

ou P t = gQtT|)m t (t l ) 
ID t (Tj) m(Ti) 
avec m(x) tel que M = [m(ij) ... m(T M )] 

sinc(-T F e B) 
et ?tt(t)= : 

_sinc((L-tF e )B)_ 



[0030] L'estimation des retards peut 6galement §tre obtenue par correlation. En considerant que le signal recu sur 
un capteur de M sources est de la forme : 



20 



25 



30 



x(t)= £a(e m )s m (t) + b(t) 

m-1 



et connaissant les angles 6 m , il est possible de considerer un signal 

s[k] = A # x[k] 

qui est forme de M versions du signal emis retardees tel que : 



(25) 



35 







s[k] = 


s(k-x 2 F e ) 








s(k-T M F e ) 



(26) 



[0031] L'autocorr elation d'une des composantes de s [k] permet d'obtenir une reference temporelle et l'intercorr6la- 
tion de cette composante avec les autres composantes donne une estimation des retards relatifs. 
as [0032] De facon similaire a la m§thode prececlente, il est aussi possible d'estimer les parametres spatiaux e m con- 
naissant les retards T m . Ceux-ci sont obtenus en calculant I'equation matricielle : 



50 



(27) 



[0033] Chacune des methodes utilisees dans ce deuxieme mode de realisation du procecle selon rinvention permet 
apre* avoir estime les angles ou les retards de trouver les autres inconnues par M recherches en une dimension. Elles 
55 ont pour avantage de repousser certaines ambigurtes, mais la seconde recherche apparatt tres dependante de la qua- 
lite de la premiere estimation. Une alternative serait d'effectuer, comme cela a 6t6 d6crit par M. J. GROUFFAUD dans 
son m6moire de these de Doctorat de TEcole Normale Superieure de Cachan de juin 1997 intitule < ( Identification 
spatio-temporelle de canaux de propagation a trajets multiples > > . des recherches successives en alternant les m§tho- 



7 
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des pour affiner la precision de Testimatton des parametres. 

[0034] Enfin, suivant un troisieme mode de realisation du proc6de selon ('invention, I'estimation des param^tres e m 
et x m a lieu de iagon conjointe en appliquant des m&hodes de determination des r6ponses impulsionnelles fondees sur 
une m d6composition en sous-espaces, un sous-espace signal de dimension M egale le nombre de sources et un sous- 
espace bruit, complementaire du sous-espace signal de dimension N(L + 1)-M, ou en appliquant la m§thode du maxi- 
mum de vraisemblance. 

[0035] Pour estimer conjointement les angles et les retards par la m£thode des sous-espaces, un projecteur P Rh est 
cre£ a partir de la matrice d'autocorrelation des filtres 



10 



15 



20 



[0036] L'estimation h du vecteur h = [h ,[0]. h ^1] h,[L] h N [0], h N [1], .... h N [L]] 1 a lieu en recherchant les 

maxima de la quantite 

P<h) Tne,) P R ,g(e,) <28) 



ou g(e,x) = a(9)«>m(T) et 



25 



M 

m=1 



avec 



30 n — rr P> ] $m 

Mm - a m e 

[0037] Dans la recherche des paramfetres par la m&hode du maximum de vraisemblance les parametres cherch£s 
35 sont les minimum du critere, 

J(ei.ti e m ,x m )=tr[l N(L+1) -G(GtG)~ 1 GtR h J 

OU G^eLT^ — ^em.Tjn)] 



40 



45 [0038] Cette recherche a lieu dans I'espace de dimension M x 2 des parametres (6 1( t 1( t m ). 

[0039] L'optimisation des differents modes de realisation du proc6d6 selon I'invention peut avoir lieu en appliquant 
une methode connue d'optimisation du type Gauss-Newton. Apres avoir trouv£ un point d'initialisation par les m6thodes 
d6crites preakientes, une optimisation des angles e m peut avoir lieu en deroulant I'algorithme suivant : 



so Initialisation : 6 0 = e |njl 
Iteration k : 

Calcul du gradient v k _., et du hessien H k _., au point 6 k ^ 
Determiner 



55 



e k = e k _ 1 -Hj: 1 _ 1 v k _., 
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Test de la convergence 

si nv^ll <c STOP sinon k= k + 1 

[0040] Le meme type d'optimisation peut etre utilise pour optimiser la recherche des retards x m . 
[0041] Un autre mode d'optimisation peut etre mis en oeuvre par une version approch6e du maximum de vraisem- 
blable suivant par exemple la methode decrite dans IEEE transaction on signal processing, Vol. 39, N1 1, Novembre 
1991 des auteurs M. Viberg, B. Ottersten, T. Kailath ayant pour titre : < < Detection and estimation in sensor arrays using 
weighted subspace fitting ) ). 



9 
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ANNEXE A 
GLOSSAIRE 

A.1. Conventions standards : 

x* : conjuguee de x 

x* : transposee du vecteur x 

x f : transposee conjuguee du vecteur x 

X 1 : transposee de ia matrice X 

X f : transposee conjuguee de la matrice X 

X" 1 : inverse de X 

X # : pseudo-inverse de X : X # =(x r x) \ r 

tr(X) : trace de X 

E[.] : esperance mathematique 

X : valeurs estimees 
x* y : convolution des signaux x et y 
5(t) : distribution de Dirac 
sinc() : fonction sinuscardinal 
® : produit de Kronecker 

A.2. Constantes entieres : 

N :nombre de capteurs composant I'antenne 

M ; nombre de trajets 

L : longueur de la reponse impulsionnelle 

K : nombre d'echantillons temporels pour composer Tobservation 
spatio-temporelle. 

A.3. Variables scalaires : 

x n (.) :observation sur le n ,eme capteur 
s() ; signal emis par la source 
b() : bruit 

a 2 : variance du bruit 

h n () : reponse impulsionnelle du n ,emo canal 
F. : frequence d'echantillonnage 
T e : periode d'echantillonnage 



10 
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am : amplitude du m' Am * trajet 

4> m : dephasage du m'* mo trajet 

6 m : azimut du m i6rne trajet 

Ym : site du m i6me trajet 

x m : retard de groupe du m ieme trajet 

a n (6) : gain du n' eme capteur dans la direction (0) 

J() : critere a optimiser 

A.4. Variables vectorielles et matricielles : 

x :observation spatiale sur I'antenne, de longueur N 

y observation spatio-temporelle sur I'antenne, de longueur NK 

h : vecteur des reponses impulsionnelles, de longueur N(L+1) 

a(6) : vecteur directionnel spatial, de longueur N 

m(r) : vecteur temporel, de longueur L+1 

l q : matrice identite de taille qxq 

R x : matrice de covariance du vecteur x 

R y : matrice de covariance du vecteur y 

Rh : matrice de covariance du vecteur h 

M : matrice des M vecteurs ternporels : M = [m(xi), m("CM)] 

A : matrice des M vecteurs directionnels : A = [a(8i) ( - a(0|^)] 

CI : matrice des attenuations des M trajets 

H : matrice des filtres mis en cotonnes. de taiile (L+1)xN 

G : matrice forrnee des reponses impulsionnelles 
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ANNEXE B 

LE MAXIMUM DE VRAISEMBLANCE DETERMINATE 



En considerant un signal de la forme : 
y k = A(6)s k + b k 

et le bruit comme etant gaussien, centre de variance a 2 , la 
distribution de probability du bruit est de la forme : 



La distribution de probability des donnees qui est a maximiser est 
donnee par la relation : 

1 — ^K-As k | 2 

f(Xk|Sk) = >r^ exp ° 



La vraisemblance des donnees pour le vecteur des donnees 
concatenees est alors : 

l(xi.x 2 x k )=n-rr^ exp ° 



k-i*'V< 



Pour maximiser la vraisemblance cela est equivalent a minimiser 



la relation 



K 

-log(L(x 1f x 2t .... x k )) = -MKIoga 2 + ]T|x k - As k | 

k=1 

La minimisation de ce critere par rapport a 0 amene a minimiser 

K 2 

e K = 2>k " As k! 
k = 1 

=>s k ^A'a) A T x k 
d'ou : 

K 

e K =Z 
k=1 



l-A^A) 1 A f jx k 
K| 2 =" k "k = t<-(uju k )=tr(u k uj) 



K 



= ZKI 

k = 1 



comme 
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e k devient : 




Ftevendications 

1. Procede de localisation spatio-temporelle d'un radio-mobile relativement a une station de base ernettrice/recep- 
trice caracterise en ce qu'il consiste : 

a calculer (1) la reponse impulsionnelle du canal de propagation entre le radio-mobile et la station emet- 
trice/receptrice, 

a estimer (2) a partir de la reponse impulsionnelle les couples (angles, retard) du trajet de propagation direct 
et des multitrajets resultats de reflexions sur des obstacles situes dans le champ de propagation des ondes 
entre la station de base et le radio-mobile, 
a calculer (2) pour chaque trajet les attenuations ctj, 

et a calculer (3) la distance d separant la station emettrice/receptrice du radio-mobile en considerant seule- 
ment les retards et les attenuations du trajet direct et d un multitrajetquelconque. 

2. Procede selon la revendication 1 caracterise en ce qu'il consiste a determiner la distance d par la relation : 



d = 




ou t 2 - est I'ecart des temps de propagation entre le trajet direct et le trajet reflechi, 
c est la celerite des ondes transmises sur le canal de propagation, 



^1 
a 2 



est le rapport des attenuations des propagations sur le trajet direct et le second trajet, 
et n la parti e due a la reflexion de I'onde sur le second trajet. 



3. Procede selon la revendication 2 caracterise en ce qu'il consiste pour calculer (2) les attenuations et og des 
ondes sur le trajet a calculer une matrice diagonale Q. def inie par le produit matriciel : 



M* HA* 



ou 




est la matrice des coefficients de la reponse impulsionnelle de canal, 
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H = [h 1 ,h 2 ,-,h N ] 

et A 1 est la matrice transposee des vecteurs directionnels a^ (e n , y n ) correspondant aux multitrajets. 

Procede selon Tune quelconque des revendications 1 a 3. caracterise en ce qu'il consists a determiner la reponse 
impulsionnelle du canal (1) au moyen d'une antenne multi-capteurs coupe a un f iltre causal, a reponse impuls.on- 
nelle finie. 

Procede selon lune quelconque des revendications 1 a 4. caracterise en ce qu'il consiste a etfectuer une estima- 
tion des angles et des retards en estimant les parametres temporels a partir de la connaissance des angles e m . 

Procede selon Tune quelconque des revendications 1 a 4, caracterise en ce qu'il consiste a etfectuer une estima- 
tion des angles et des retards en estimant les parametres spatiaux 6 m a partir des retards T m . 

Procede selon Tune quelconque des revendications 1 a 4, caracterise en ce qu'il consiste a etfectuer une estima- 
tion conjointe des angles 0 m et des retards x m . 

Procede selon I'une quelconque des revendications 1 a 7, caracterise en ce qu'il consiste a estimer les 
parametres temporels et les directions d'arrivee des multitrajets par la mise en oeuvre de I'algonthme connu sous 
le nom de < ( MUSIC )). 

Procede selon Tune quelconque des revendications 1 a 7. caracterise en ce qu'il consiste a estimer les parametres 
temporels et les directions d'arrivee des multitrajets par la mise en oeuvre de Talgorithme determinate connu sous 
la denomination de ( < maximum de vraisemblance )). 
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